Ratios of Multijet Cross Sections in pp Collisions at y/s = 1.8 TeV 
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We report on a study of the ratio of inclusive three-jet to inclusive two-jet production cross 
sections as a function of total transverse energy in pp collisions at a center-of-mass energy y^ =1.8 
TeV, using data collected with the D0 detector during the 1992-1993 run of the Fermilab Tevatron 
Collider. The measurements are used to deduce preferred renormalization scales in perturbative 
C(c*s) QCD calculations in modeling soft-jet emission. 



A primary manifestation of Quantum Chromodynam- 
ics (QCD) in pp collisions at a high center-of-mass energy 
[s/s = 1.8 TcV) is the production of jets with large trans- 
verse momenta. Typically, the hard interaction of par- 
ton constituents of a proton and an antiproton produce 
two hard back-to-back jets. However, a fraction of the 
time, additional jets are also produced. In the absence of 
an all-orders QCD calculation, jet production rates as a 
function of jet energy are predicted by fixed-order calcu- 
lations in perturbative QCD (pQCD). In this paper, we 
investigate the dependence of these calculations on the 
choice of parton distribution functions (pdf) and partic- 
ularly renormalization and factorization scales 

We examine the ratio of inclusive three-jet production 
to inclusive two-jet production, which reflects the rate 
of gluon emission in QCD jet production processes. A 
three jet cross section explicitly offers the opportunity 
to investigate a scale difference at a secondary vertex. 
Taking the ratio reduces systematic uncertainties. 

Although this issue has inherent theoretical interest, it 
is also important because QCD multijet production is fre- 
quently a background to rare processes: phenomenolog- 
ically confirmed prescriptions for renormalization scales 
are essential for predicting background rates and for de- 
signing efficient triggering schemes for rare processes at 
future colliders m. Lastly, when higher order QCD cal- 
culations become available, this ratio may be useful for 
providing another accurate measure of the strong cou- 
pling constant a^. 

The data used in this analysis, corresponding to an 
integrated luminosity of sa 10 pb^ , were recorded dur- 
ing the 1992-1993 Tevatron collider run. The D0 de- 
tector is described in detail elsewhere 0. Jet detection 
primarily utilizes the uranium-liquid argon calorimeters, 
which have full coverage for pseudorapidity jryl < 4 where 
rj = — ln[tan(0/2)] and is the polar angle relative to 
the direction of the proton beam. Initial event selection 
occurred in two hardware trigger stages and a software 
stage. The first hardware trigger selected an inelastic pp 
collision as indicated by signals from trigger hodoscopes 
located near the beams on either side of the interaction 
region. The next stage required transverse energy above 
a preset threshold in calorimeter towers of 0.2 x 0.2 in 
A77 X A0, where is the azimuthal angle. Selected events 
were digitized and sent to an array of processors. Jet can- 
didates were then reconstructed with a cone algorithm 
and the event recorded if any jet transverse energy {Et) 
exceeded a specified threshold. Five such inclusive trig- 
gers had thresholds of 20, 30, 50, 85, and 115 GeV. 

Jets were reconstructed offline using an iterative fixed- 
cone algorithm with a cone radius TZ — 0.7 in rj ~ (j) 
space. The Et of each jet was corrected for effects due to 
the underlying event, additional interactions, noise from 
uranium decay, the fraction of particle energy deposited 
outside of the reconstruction cone, detector uniformity, 
and detector hadronic response. A discussion of the jet 
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FIG. 1. The ratio R32 as a function of Ht for Et thresh- 
olds of 20, 30, and 40 GeV (|?7jot| < 3). Error bars indicate sta- 
tistical and uncorrelated systematic uncertainties, while the 
distribution at the bottom shows the correlated systematic 
uncertainty for the 20 GeV threshold. 



algorithm, energy scale calibration and resolution can be 
found inRefs. |-||. 

We measure the ratio of the inclusive three-jet to the 
inclusive two-jet cross section 
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(T3 cr{pp -^ n jets + X ; n > 3) 
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as a function of the scalar sum of jet transverse energies 
{Ht — X] ^T ) ■ ^^^ measurement is performed for four 
distinct sets of selection criteria for all jets in the event: 
Et thresholds of 20,30, or 40 GeV for |r?jet| < 3, and 
Et > 20 GeV for |77jet| < 2. Three thresholds were cho- 
sen to study threshold dependence, and the minimum 
threshold was chosen to maximize statistics for which 
jet reconstruction efficiency was nearly 100%. Both in 
the data analysis and in the QCD calculation, a jet con- 
tributes to Ht and to the jet multiplicity if it passes all 
selection criteria and satisfies the Et and r/jct require- 
ments. 

Figure |l| shows the ratio R32 as a function of Ht for 
Et thresholds of 20, 30, and 40 GeV for JTyjotl < 3. The 
five trigger samples listed in the figure contribute in sep- 
arate regions of Ht, as indicated by the symbols. The 
distribution at the bottom of the figure shows the corre- 
lated systematic uncertainties for the 20 GeV threshold. 
This uncertainty is the maximum offset in the ratio ob- 
tained by a one standard deviation change in the correc- 
tion to the jet energy scale. Error bars indicate statistical 
uncertainties (calculated using the appropriate binomial 
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TABLE I. Values of R:j2 with their uncorrelated and correlated uncertainties for the indicated jet Et threshold and r/jet 
criteria. Uncorrelated uncertainties include statistical and uncorrelated systematic uncertainties added in quadrature. 



prescription for a statistically correlated ratio) as well 
as uncorrelated systematic uncertainties arising from all 
selection criteria. Table ffl displays the measurements in 
bins of Ht, showing the uncorrelated and correlated un- 
certainties for the four selection criteria. 

JETRAD is a next-to-leading-order Monte Carlo gen- 
erator for describing inclusive multijet production. The 
generated 2-jet and 3-jet events are inclusive, and there- 
fore the ratio of these cross sections should be equivalent 
to the measured i?32. CTEQ4M Q pdf are used in the 
JETRAD simulations. The jet finding algorithm in je- 
TRAD approximates the algorithm used in D0 data re- 
construction. Jets generated by jetrad are individually 
smeared according to known detector resolutions. Two 
partons are combined if they are within 7?.scp — 1.37?., as 
motivated by the separation of jets in the data [B| and, 
just as in the data, a jet is included if its Et and 77jct 
meet the chosen selection criteria. 

In pQCD, the renormalization procedure introduces a 
mass scale fin to control ultraviolet divergences in the 
calculations. A factorization scale fip, introduced to han- 



dle infrared divergences, is assumed to be equal to fiR in 
all predictions described in this paper. QCD provides the 
evolution of as with fin, but not its absolute scale. Unless 
otherwise indicated, the renormalization scale fiR = XHt 
will be used for the production of the two leading jets, 
where the constant A, the coefficient of the hard scale, 
will have a nominal value of 0.3, but will be allowed to 
vary as described below. To study the possibility of hav- 
ing a different scale for the production of additional jets, 
the renormalization scale of the third jet is varied from 

/i]^ = XHt (same as for the leading jets) to a scale pro- 
fs) (3) 
portional to the Et of the third jet fi)^ ex E"^ . Also, a 

scale proportional to the maximum jet transverse energy 
l^^max-j jg studied, as this is a standard form used for 
comparisons of jetrad to measured jet cross sections. 

Figure || shows the measured R32 as a function of Ht 
for jet Et > 20 GeV and JTyjotl < 2. The 20 GeV threshold 
has good sensitivity to scale in the jetrad prediction and 
has reduced statistical uncertainty. The central rapidity 
region has the best understood jet energy uncertainties 
and correlations. The plot contains four smoothed distri- 



cr 



0.8 F 

■O 

'0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




- Et>20 GeV 




<2 




^-.-i,^.- 




.2__^ 1 1 


■ MeasL 


"T — """^^^^^^-^^^ - ■ 1 


jred R32 


: ;';# 




JETRAD 


/i,"'^' = 0.3 H, 


'-^r 




— Mr = 


0.3 Ht 


y- 






= 1.0 ^r 

= 2.0 e;^' 






JETRAD 








-.-.- ^R = 


0.6 Er™ 


: Data correlated 
-^-____^_^ uncert 


ainty 







100 



200 



300 



400 500 600 

Ht (GeV) 

FIG. 2. The ratio i?32 as a function of Ht, requiring 
jet Et > 20 GeV and |ryjot| < 2. Error bars indicate sta- 
tistical and uncorrelated systematic uncertainties, while the 
histogram at the bottom shows the correlated systematic un- 
certainty. The four smoothed distributions show the jetrad 
prediction for the renormalization scales indicated in the leg- 
end. 



butions corresponding to jetrad predictions for the fol- 
lowing renormalization prescriptions (shown for A = 0.3): 

• /ifl — XHt for the two leading jets, 

— fiji — XHt also for the third jet (solid) 

(3) _ R.(3) 



,(3) 



E}p ' for the third jet (dashed) 



r(3) 



Mi? 



M)^^ = 2E)j:' for the third jet (dotted) 
: 0.6£:j'^'^ for all jets (dash-dot). 



All predictions demonstrate the same qualitative behav- 
ior as the i?32 measurement, that is, a rapid rise below 
Ht = 200 GeV (associated with the kinematic thresh- 
old), a leveling off, then a slight drop at highest Ht (as- 
sociated with the reduced phase space for additional ra- 
diation for high Et jets). Although jetrad predictions 
for the ratio are found to be insensitive to the choice of 
pdf, they do depend on the choice of TZsep- Allowing 7?.scp 
to vary such that neighboring jets are all merged or all 
split causes a 3% decrease or increase in the ratio, re- 
spectively, with only a slight effect on the shape of the 
distribution in Ht- 

For a quantitative comparison, we use a x^ covariance 
technique, defining 

n n 
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FIG. 3. x^/dof as a function of A, comparing data to je- 
trad predictions for several renormalization prescriptions for 
the Et > 20 GeV, \rijct\ < 2 selection criteria. 



where Di and Ti represent the i'^ data and theory el- 
ement, respectively, and C"-*- is the inverse of the co- 
variance matrix. This matrix incorporates uncorrelated 
uncertainties in the measurement and statistical uncer- 
tainties in the simulation, with correlated uncertainties 
included for the absolute jet energy in the data and for 
the uncertainty from resolution smearing in jetrad (not 
shown explicitly in Fig. g). Although some of the pre- 
dictions do not visually overlap with the data, accept- 
able agreement is found for some scales because of the 
strong point-to-point correlations of the data uncertain- 
ties which are taken into account in the x^- Figure y 
shows the x^ per degree of freedom (x^/dof) as a function 
of the parameter A, for the Et > 20 GeV, |77jct| < 2 selec- 
tion criteria. The degrees-of-freedom equal the number of 
data points (28). The horizontal line indicates the x'^/dof 
obtained using the A independent scale ^ir = 0.6Et^^ for 
all jets. This scale yields good agreement with measure- 
ment (probability p > 57%) for the Et > 20 GeV crite- 
ria, but the x^ rises (and the corresponding probabilities 
decrease) for the higher Et thresholds (not shown). 

For A-dependent scales, the best fit is specified by the 
A that minimizes the x^ ■ The scales proportional to Ej. 
for the third jet do not provide a good fit {p < 5%) for 
any A, as seen in Fig. |3|. While there is fair agreement 
in the wider region of pseudorapidity jryjotl < 3 for cer- 
tain regions of A (not shown), these do not correspond 
to the same values for different Et thresholds, making 
the applicability of this scale prescription unsuitable for 
predicting production rates for additional jets. 

The JETRAD prediction assuming a scale iir — XHt 



provides the best description of the data for A between [8] B. Abbott et al, "Fixed cone jet definitions in D0 and 
0.30 and 0.35 (p > 80%). Moreover, the x^ is also mini- R(sep)", Fermilab-Pub-97-242-E, unpublished, 

mized in the A « 0.30 region for the other selection crite- 
ria (not shown) making this scale choice the most robust 
of all the /ifl scales studied. 

In conclusion, wc have measured the ratio of the in- 
clusive three-jet to the inclusive two-jet cross section as 
a function of total scalar transverse energy Ht and com- 
pared the results to jetrad predictions. The greatest 
sensitivity to the choice of renormalization scale is for 
the lowest Et threshold of 20 GeV. Although no predic- 
tion accurately describes the ratio through the kinematic 
threshold region, a single iir scale assumption in the cal- 
culation for all jets is found to adequately describe the 
rate of additional jet emission when correlated uncertain- 
ties are accounted for in a x^ comparison. Specifically, 
a scale of hr = XHt for all jets, where A = 0.3, yields 
a prediction consistent with the measurement for all jet- 
selection criteria examined. A scale of ^ji = O.QEi^'^^ for 
all jets also provides a sufficient description at the lowest 
jet Et threshold. The introduction of additional scales 
does not significantly improve agreement with the data. 
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